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Abstract. We use a hydrodynamical model to describe the evolution of the collision system at collision energies 
■s/s =130 and 200 GeV. At lower -y/s =130 GeV energy we compare the results obtained assuming fast or slow 
thermalization (thermalization times To = 0.6 and 4.1 fm/c, respectively) and show that slow thermalization fails 
to reproduce the observed anisotropy of particle distribution. At ^/s = 200 GeV collision energy our results show 
anisotropies similar to those observed at yfs =130 GeV. 

INTRODUCTION 

The azimuthal anisotropy of particle distributions in heavy ion collisions is a signal of rescattering among the particles. 
The larger the anisotropy, the stronger the rescattering. Since the hydrodynamical models assume zero mean free path 
and thus infinite scattering rate, they provide an upper limit to observable anisotropies. In our earlier papers [1, 2, 3] 
we have argued that since the observed anisotropies are well reproduced by a hydrodynamical model, the collision 
system must thermalize fast. Here we constrain our model further by requiring it to reproduce the observed single 
particle /^-distributions [4] and study the consequences of such a long thermalization time that there is no thermal 
plasma but the hydrodynamical expansion begins in hadronic phase. We also tune our model to reproduce the charged 
particle multiplicity at ^fs = 200 GeV collision energy published recently by the PHOBOS collaboration [5] and show 
our results for the charged particle anisotropy and pion and antiproton /vspectra at the higher energy. 

HYDRODYNAMICAL MODEL AND ITS INITIALIZATION 

To reduce the complexity of the numerical calculations we assume boost invariant scaling flow in our hydrodynamical 
model and solve the equations numerically only in the two transverse dimensions. We also apply ideal fluid hydrody- 
namics and assume that viscous effects are negligible. Depending on the initialization we apply two different equations 
of state: one with a phase transition at T c = 165 MeV or purely hadronic equation of state, called EoS A and H in [1], 
respectively. 

In our model we assume that the chemical and kinetic freeze-out take place at the same temperature. In other words 
the system maintains both kinetic and chemical local equilibrium during the entire evolution from thermalization to 
freeze-out. This makes it impossible for us to fit both proton and antiproton yields simultaneously since the thermal 
model studies [6] point to chemical freeze-out temperatures much higher than the kinetic freeze-out temperature 
Tf = 100 - 140 MeV required to produce sufficiently flat /?,-spectra. We counter this problem by the following 
approximation: We fix our initial baryon density to reproduce the observed p/p ~ 0.6 [7] ratio at T — 165 MeV, 
i.e. immediately after the system has hadronized. We then let the system evolve to its freeze-out temperature and scale 
the calculated proton and antiproton yields to their values at T = 165 MeV by hand. This scaling factor depends on 
the kinetic freeze-out temperature and has to be adjusted accordingly. 



Fast thermalization 



We have discussed different ways to parametrize the initial state of hydrodynamic evolution in [3]. It was found that 
neither a parametrization where the density (energy or entropy) is proportional to the local number of participants or 



binary collisions is sufficient to reproduce the observed centrality dependence of multiplicity [8]. A combination of 
these like the two component model presented in [9] is required instead. 

To obtain the results presented here we have used a combination of parametrizations eWN and eBC described 
in [3] 1 . In this case we assume the initial energy density to be proportional to the scaled sum of number of participants 
and binary collisions per unit area: 

e(s,T ;b) ^ xK e (x )n^ N (s;b) + (I -x)£ e (x )n B c(s;b), (1) 

where «wN(s;b) is the number of participants and «Bc(s;b) the number of binary collisions per unit area in the 
transverse plane. The normalization constants K e (io) and K e (xo) are the same than in [3] and the mixing parameter 
x = 0.4 was found to reproduce the observed centrality dependence of multiplicity. The equation of state with phase 
transition (EoS A) has been used to calculate results with short thermalization time (Xq = 0.6 fm). 



Slow thermalization and free streaming 

The parametrization described above is reasonable if thermalization time is sufficiently small and produced particles 
do not have time to move far away from the place where they were produced. On the other hand, if it takes several 
fermi for the system to thermalize, the particles will move during thermalization and smear any structure in the system. 

To estimate what the shape and distributions of the system might be after a long thermalization process, we 
developed further the ideas outlined in [ 1 0] . We approximate the expansion of the system during thermalization by free 
streaming. In this approach the number density of initially produced excitations is again taken to be proportional to 
the scaled sum of number of participants and binary collisions whereas the momentum distribution is taken to be that 
of pions at temperature of T = 200 MeV. The free streaming is allowed to continue until thermalization time To when 
energy density distribution due to free streaming particles is calculated and taken to be the thermal energy distribution. 
In this approximation we do not assign any radial velocity to the system at the end of free streaming. The effects of 
build-up of flow during thermalization will be discussed in [11]. 

The distribution of particles in a free streaming system is given by [12] 

«(x,P,T ) = fo(jft,p' z )S±(s-y_ L (z -'Zi),b), (2) 

where 5j_(s,b) °< xK„n^(s;b) + (1 — x)/T/,«Bc(s;b) gives the density distribution of the initial excitations and 
fo(Pt,Pz) i s Boltzmann distribution. Of the parameters p' z = XqP z /%i is the scaled longitudinal momentum which 
takes longitudinal free streaming in boost invariant system into account, To is thermalization time and x, is the time 
when the distributions before free streaming are specified. We use T; = 0.6 fm and fix thermalization time by requiring 
that maximum temperature in central collisions is low enough for hadrons to exist. We choose this somewhat arbitrary 
limit to be T = 200 MeV. This corresponds to e — 2.3 GeV/fm 3 in the hadronic equation of state (EoS H) we use 
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FIGURE 1. Initial energy density distribution in b — 2.4 fm collision scaled by maximum energy density (left) and initial spatial 
deformation of the system as function of impart parameter (right) using parametrizations for fast and slow thermalization. 



1 A combination of sWN and sBC of [3] is equally possible. 
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FIGURE 2. p ( distributions of negative pions (left) and antiprotons (right) at ^/s =130 collision energy compared with pre- 
liminary PHENIX data [4], The data and curves correspond to centralities a/a tot = 0-5%, 5-15%, 15-30%, 30-60%, 60-92% and 
impact parameters b = 2.4, 5.1, 7.0, 9.6 and 12.1 fm, respectively (top to bottom). These data sets are successively scaled down by 
10"",n = 0,1,2,3,4. 



with this initial state. After the remaining parameters K n , Kj, and x are chosen to reproduce the observed centrality 
dependence of multiplicity [8], this procedure leads to a thermalization time To = 4.1 fm. 

As shown in Fig. 1 this parametrization leads to a larger system which density gradients are flatter and the spatial 
deformation smaller than in the fast thermalization parametrization (Eq.(l)). 



COMPARISON TO OBSERVABLES 
p t - distributions 

After fixing the proportionality constants K e and K e and mixing x in Eq.(l) the only remaining parameter in our 
model is the freeze-out temperature. We fix it to reproduce the preliminary /^-distributions of pions and antiprotons in 
most central collisions measured by the PHENIX Collaboration [4]. This requirement leads to very different tempera- 
tures. If early thermalization is assumed, 7/ w 140 MeV leads to a nice fit. On the other hand long thermalization time 
and free streaming of particles during thermalization smoothens the pressure gradients so much that a low freeze-out 
temperature of Tf w 100 MeV is needed to create large enough transverse flow (Fig. (2)). This difference in freeze-out 
temperatures leads also to very different scaling factors for the antiproton yield, 2 and 27 for To = 0.6 fm and To = 4.1 
fm, respectively. After scaling the yields to their values at T c h = 165 MeV, we reach an acceptable reproduction of 
measured spectra in most central collisions. 

Having the parameters of our models fixed we can now compare our results to spectra of non-central collisions [4] 
(see Fig. 2). The average number of participants in five different centrality classes [13] set our choice of impact 
parameter in these bins to 2.4, 5.1, 7.0, 9.6 and 12.1 fm, respectively. The early thermalization assumption leads to 
very good agreement with the preliminary data even for the most peripheral collisions. On the other hand, the late 
thermalization approach performs adequately in central and semi-central collisions, but leads to far too steep slopes 
in peripheral collisions. This failure could be expected since if one assumes a long thermalization time even in central 
collisions, it is probable that the system does not thermalize at all in peripheral collisions and our model fails. 



Azimuthal anisotropy 

So far we have used different observables to fix parameters of our model and the real test is in the comparison 
to the observed anisotropics of particle distribution [14, 15]. As we have reported in our earlier papers [1, 2, 3], 
hydrodynamical model leads to an excellent fit to observed data if thermalization time is short. We show our present 
results in Fig. (3). The present parametrization (To = 0.6 fm) leads to a slightly smaller differential anisotropy for 
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FIGURE 3. Differential anisotropy of pions and protons+antiprotons (upper and lower curve, respectively) in minimum bias 
collisions compared to preliminary data [14]. 



pions and larger anisotropy for antiprotons than shown in our previous works. Both changes are mainly due to higher 
freeze-out temperature (see [2]). The deviation from data is still within systematical error and as reported in this 
conference [16], more careful analysis of the data to make systematic errors smaller may lead to slightly smaller V2- 
Also we want to remind that our present fit to freeze-out temperature Tf f« 140 MeV is preliminary and may change 
when we fit more data and explore different parameter combinations. 

On the other hand the results for long thermalization time approach are well below the data and show that - at least 
in the present formulation - the free streaming period changes the shape of the initial system sufficiently to prevent 
the build-up of observed anisotropies. This corroborates our earlier argument that early thermalization is necessary to 
achieve the observed large anisotropies in particle distributions. 



COLLISIONS AT Vs = 200 GEV 

We use the charged particle multiplicity dN/dr\ = 650 ± 35 in 6% most central collisions measured recently by the 
PHOBOS Collaboration [5] to fix our parametrization for yfs = 200 GeV collisions. We assume that all the other 
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FIGURE 4. p t distributions of pions and antiprotons at 5% most central ~Js = 130 and 200 GeV collisions (left) and differential 
anisotropy of charged particles in minimum bias sfs = 130 and 200 GeV collisions compared with preliminary STAR data at 
sfs = 130 GeV [15] (right). 



parameters like freeze-out temperature (Tf w 140 MeV) and thermalization time (To = 0.6 fm) stay unchanged. We do 
not have any method to calculate the change in baryon stopping and thus the change in initial baryon density. Therefore 
we ignore the question of the change in stopping altogether and use the same initial baryon density at both energies. 
For the same reason we do not scale the calculated antiproton yield at freeze-out either but show the result as it is to 
demonstrate the change in slopes. 

The calculated pion and antiproton distributions at \/s = 130 and 200 GeV are shown in Fig. (4). The changes in 
slopes are small but still large enough to increase E t by ~ 18% even if the observed increase in particle multiplicity is 
14% [5]. Our calculation leads to the increase of the average energy density at x = 1 fm from 5.4 to 6.8 GeV/fm 3 and 
the increase of the maximum temperature (at Xo = 0.6 fm) from 355 to 375 MeV. If the changes in /?,-spectra are small 
they are even tinier in anisotropies. In Fig. (4) differential anisotropy for charged particles in minimum bias collisions 
is shown with preliminary STAR data for y/s = 130 GeV collision. The change is negligible and it will be interesting 
to see whether the data deviates from the calculation around p t ~2 GeV also at \/s = 200 GeV collision energy. 



SUMMARY 

We have shown that it is possible to fit both the experimental p t spectra and anisotropy of particle distribution if 
short thermalization time and hydrodynamical behaviour of the collision system is assumed. On the other hand long 
thermalization time allows the system shape change too much before the build-up of flow begins and generates 
sufficient anisotropies. This corroborates our earlier argument that fast thermalization is required to explain the 
experimental data. 

We have also made the first calculations at y/s = 200 GeV collision energy. Our results show only slightly flatter 
spectra and basically no change in differential V2 when compared to results at lower y/s = 130 GeV energy. 
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